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Abstract 

Using a first principles based, magnetic tight-binding total energy model, the 
magnetization energy and moments are computed for various ordered spin configu- 
rations in the high pressure polymorphs of iron (fee, or 7-Fe, and hep, or e-Fe), as 
well ferromagnetic bcc iron (a-Fe). For hep, a non-collinear, antiferromagnetic, spin 
configuration that minimizes unfavorable ferromagnetic nearest neighbor ordering 
is the lowest energy state and is more stable than non-magnetic e iron up to about 
75 GPa. Accounting for non-collinear magnetism yields better agreement with the 
experimental equation of state, in contrast to the non-magnetic equation of state, 
which is in poor agreement with experiment below 50 GPa. 
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1 Introduction 



Magnetism is known to be important in the phase stability, structure and 
elastic properties of iron. For example, a-Fe, the ground state at ambient con- 
ditions, would be mechanically unstable if it were not magnetic. Even above 
the Curie temperature, Tc there are local magnetic moments in a-Fe. Face- 
centered cubic iron (fee or 7-Fe) has incommensurate magnetic correlations 
which change rapidly with volume, arid give rise to the anti-Invar effect (large 
thermal expansivity) (jMrvasov et al. . 1992f ). The magnetic behavior of hep 



iron is important for high pressure materials research, and for interpreting 
high pressure experiments aimed at understanding the Earth's inner core. Hep 
iron is not quenchable to zero pressure, so magnetic studies must be made in 
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situ at high pressures. e-Fe was long thought to be non-magnetic, due to sev- 
eral Mossbauer experiments that showed no hyperfine splitting in e-Fe, even 
down to helium temperatures. However, self-consistent first-principles compu- 
tations s how a magnetic ground stare for e-Fe which is stable up to about 
50 GPa ( Steinle- Neumann et al. . 1999[ ). Here non-coUinear magnetism in e- 



Fe is explored using a magnetic tight-binding model fitted to first-principles 
calculations. 

When magnetic moments are coUinear, electrons can be considered to be "spin- 
up" or "spin-down" in a global sense, i.e. there is a global magnetic quantiza- 
tion direction. This means that one can solve for the spin-up and spin-down 
electrons separately, and then combine the results to compute the total charge 
density. When moments are oblique to each other, the spin state is said to be 
"non-collinear." In that case the problem does not factorize, and one must 
diagonalize a Hamiltonian of twice the order of the coUinear case. A system 
can be non-collinear either in response to chemical or thermal disorder, or in 
order to minimize frustration. Both hep and fee lattices are frustrated with 
respect to antiferromagnetism, in that one cannot order these lattices anti- 
ferromagnetically so that all neighbors have opposite spins. This is known to 
lead to spin-waves and non-collinear magnetism in fcc-Fe. Also, a common ex- 
ample of non-collinear spins occurs on heating bcc-Fe above the Curie point, 
where the spins disorder dynamically. A comprehensive review of the theory 
of non-collinear magnetism is given in ISandratskii 



2 Method 



To study non-colline ar magnetism in Fe , a first-principles based non-m agnetic 



tight-binding model (|Cohen et al.l . ll994ll997l:IWasserman et al.l.ll99fit) 



IS com- 



bined with a model for magnetism ( Pickettl . 19961 : iMukheriee and Coh^ , 2001 ) 
The Hamiltonian (H) of the system is given by. 



H = Hn -\- R H, 



coUinear 



rR, 



where Hq is the (doubled) non-magnetic Hamiltonian, Hcoiiinear is the coUinear 
magnetic Hamiltonian and R is the rotation matrix for spin directions, as de- 
scribed below. The overlap matri x 5* is unchanged from the non-magnetic non- 
ortho gonal tight-binding model ( Cohen et al. . 19941 1997 : Wasserman et al. . 
1996t ). Hq was fit to eigenvalues and total energies from an extensive set of 



non-magnetic Linearized Augmented Plane Wave (LAPW) results within the 
Generalized Gradient Approximation (GGA) dPerdew et al. . 1996f ) in the bcc, 
fee, and hep str uctures as functions of pressure and strain, and ha s been ex- 
tensively tested (jCohen et al.l . ll994L Il997l : IWasserman et al.l . ll996h . 
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The magnetic Hamiltonian is twice the size of the original non-magnetic 
Hamiltonian, with the upper left block being spin-up, the lower right block is 
spin down, and the off-diagonal blocks are from coupling between up and down 
spins that give rise to non-coUinear magnetic solutions. Hq is doubled from 
the original non-magnetic Hamiltonian, with identical diagonal blocks, and 
zeros on the off-diagonal blocks. The only non-zero elements in the coUinear 
magnetic Hamiltonian '^coiHnear are the diagonal onsite d elements. The lat- 
ter are given by — /mj/2 for the spin-up block, and +Imi/2 for the spin-down 
block, where rrii is the magnetic moment of the atom i, and / is the St oner pa- 
rameter, which controls the strength of the exchange splitting (jPickett, 1996). 
R is the rotation matrix which dep ends upon the a tomic spin direction with 



respect to a global reference frame (jUhl et al.l . 11994 ) . The diagonal of H^ag is 



spin-up and spin-down, and the off-diagonal bands are the coupling between 
up and down, which arise from the rotation operation. Similar models have 
been us ed previously fo r incorporating magnetism within a tight-binding ap- 
proac h (|You and Heinel . Il982l : IPicket"tl . I199(tI : IPrevss et al.1 . Il997l: iMehl et al. . 
20011 ). but the magnetic behavior of e-Fe has not been addressed, nor has a 



first-principles tight-binding model previously been developed that includes 
constraining fields. 

The full spin density matrix is computed, and the constraining fields are com- 
puted self-consistently so that spin directions as well as spin moments are 
self-consistent. The spin-density matrix is given by 



Po 



"^iaSijlpja' 



(2) 



where Sij is the overlap matrix and ijjia are the eigenvectors for atom orbital i 
and spin a in the local coordinate system on each atom. In the local coordinate 
system the spin on each atom is diagonal and can be represented as up or 
down. The spin-1/2 rotation matrix elements of R that rotates from the local 
coordinate system to the global system with oblique spins is given by 



R 




cos ttC 2 sm 2 



— sm|e 2 cos|e 2 



(3) 



where 6 is the polar angle and is the azimuthal angle. The matrix elements 
given here are used to build the full rotation matrix; the only non-zero elements 
are diagonal in atom orbital indices. Without spin-orbit coupling, there is no 
coupling to absolute directions relative to the crystalline lattice, and only the 
relative angles between spins are important. Spin-orbit interactions are not 
included here; their energetic importance is insignificant compared with the 
exchange energy changes for iron. 
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The magnetic moments are given by the spin density on each atom summed 
over the orbitals. The expressions are 



e = aTccos(^ — ^) (5) 



( m 

(f) = ^\ogp^l. (6) 

The magnetic moments are determined self-consistently. One starts with an 
initial guess, which gives the exchange sphttings, ±/mj/2, in the Hamiltonian, 
find H, and diagonahze the generahzed eigenproblem. States are occupied up 
to the Fermi level, and the output moment on each atom is found from the 
weighted eigenvectors and the expressions above. The process is repeated until 
self-consistency is obtained. 

In general the output moment direction, as well as magnitude, will differ from 
the input. In this way a self-consistent process can be used to find the moment 
magnitudes and directions that are locally stable. In some cases one wants to 
compute the energy of a given magnetic structure. In order to guarantee that 
the output magnetic structure is the desired structure, it is necessary to apply 
magnetic fields that force the output moments to form the desired structure. 
Such a procedure has been implemented using a self-consistent procedure with 
the option to constrain the moments in magnitude and direction (constraining 
fields contain longitudinal components), or to allow the magnitude to adjust 
self-consistently, but force the direction to remain as desired (constraining 
fields are transverse). The procedure used is similar to, but different from, the 



procedure outlined in lUifahissv et ID (|l999>) . 



The total energy is given by 

Etotal =E, + IJ2 "^i /4 + (7) 



where Ej, is the band energy (sum of the eigenvalues). There are no additional 
potential terms; the total energy is obtained from the band structure and mo- 
ments. The second term in Eq. 7 corrects for double counting of the exchange 
interaction (^Pickett . and the third term corrects for double counting of 

the interaction with the applied staggered fields bi. 



The model differs from the conventional Stoner model ( Stonerlll938l ). in which 
a ferromagnetic instability is predicted by the inequality IN{0) > 1, where 
A^(0) is the non-mag netic density o f states at the Fermi level, and the ex- 
tended Stoner model (|Kraskol . lT987h . in which N is replaced by the effective 
density of states N{M) = M/Se, where 6e is the exchange splitting, in that 
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our model allows for different hybridization depending on magnetic state, ac- 
counting thereby to the actual magnetic structure, i.e. whether the system 
will be ferromagnetic, antiferromagnetic, or non-collinear. 



3 Results and Discussion 



3.1 bcc 



The model is in good agreement with previous self-co nsistent calculations for 
Qi-Fe, especially when a volume dependent / is used flMukheriee and Cohen . 
2001[ ). The lowest energy state for a- Fe for the model at ambient and at higher 
pressures is ferromagnet ic, in agreemen t with the first principles LAPW self- 



consistent calculations (|Stixrude et al.l . Il994^ . Moreover, even with a single 



value of volume independent / the model calculation describes well the pres- 
sure dependence of the magnetic moment and the magnetization energy. As 
pressure is increased the magnetization energy decreases smoothly (Fig. 1). 
The value of / needed for quantitative agreement is almost constant (around 
1 eV) for larger atomic volumes or lower pressures but increases for higher 
pressures. This is consistent with t he increase in the exchange interaction at 
higher pressures noted previously ( Asada and Terakura . 19921 ). Even though 
the exchange interaction increases with pressure, the non-magnetic density of 
states near the Fermi energy decreases, giving rise to a net reduction in the 
magnetic moment and magnetization energy, and finally the loss of magnetism 
at around atomic volume of 40 bohr^. 



The Stoner parameter I was adjusted at each volume to give agreement with 
self-consistent full- potential Linearized Augm ented Plane Wave (LAPW) mag- 
netization energies ( Steinle-Neumann et al.l . ll9^ 99) for each volume for bcc. / is 
expected to increase as the electronic density increase s with decreasi ng atomic 
volume, similar to what has been seen before for 7-Fe ( Kraskol . 1987 ) and a-Fe 
(Fig. 1). Fitting the results to a polynomial valid between V =50 bohr^ and 
90 bohr^ gives I{eV) = 3.4126 - 0.08583V" + 0.00103^^ _ 4.166610-^1/3 _ 



3.2 fee 



7-Fe shows a great richness in non-collinear magnetic structures as volume is 
varied, a nd here it is not exp lored in detail; much study has been done, as re- 
viewed in ISandratskii ( 1998| ). The energies and moments for fcc-Fe with spiral 
spins along (001) and 6 = tt/2 were determined using eight-atom supercells, 
allowing calculations for (OOg) with q = 0, 7r/4, 7r/2, 37r/4,and tt. A k-point 
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Fig. 1. Magnetization (a) and magnetization energy ( b) for bcc iron as fun ctions 
of volume. The self-consistent LAPW results are from IStixrude et al.l ()l994l 'l. The 
tight-binding (TB) model results are shown for constant 1=0.95 eV, and for varying 
/, with / chosen to best fit the LAPW magnetization energies. The best-fit / is also 
shown in (a). 

mesh of 12x12x4 was used, giving 72 k-points in the irreducible wedge for 
tetragonal symmetry. For 7-Fe, the complex behavior of magnetic moments 
and the magnetization energy with increasing pressures and varying spiral 
spin density wave states is qualitatively repro duced by the rnodel (Figs. 2 and 



3) compared with self-consistent calc ulations (iMrvasov et al.l . 119921 : lUhl et al, 
1994 : Sjostedt and Nordstrom . I2nn2t ). Quantitatively the results are sensitive 



to the value of /, and to get good agreement with self-consistent calculations 
it seems a smaller I, about 0.94 eV, is required than derived from bcc-Fe 
(0.99-1.02 eV for the volumes considered for fee). In any case there is some 
variation in self-consistent results for e-Fe, due to extreme sensitivity of the 
magnetic structure to basis set, k-point samphng, etc., due to the small en- 
ergy scale. Furthermore, the ma gnetic ground state in 7-Fe is se nsitive to the 
atomic momeut approximation k,^>.^>. NorH.trh ^ EHi in the tight- 
binding model; that is the moment is really a field that varies with position 
in the crystal, and is not constant on each atom. Within the atomic moment 
approximation our results are reasonably consistent with the self-consistent re- 
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suits. More det ailed comparisons with the comp rehensive non-coUinear LAPW 
calculations of ISiostedt and NordstromI ( 20021 1 inclu ding k-point convergence 
tests are called for, but have not yet been done. In ISiostedt and Nordstrom 
fj2002h 4000 k-points were used in the full Brillouin zone, compa red to our 4608 



which seems comparable. It is more difficult to compare with iMrvasov et al 



since a real space multiple scattering approach and a muffin-tin poten- 
tial approximation was used. 



3.3 hep 



Previous collinear first principles c alculations show stability of a n antiferro- 
magnetic state, afmll, in hep iron ( Steinle- Neumann et al. . 1999| ). Moreover, 



the computed equation of state of afmll greatly improves agreement with the 
experimental equation of state. An hep lattice is frustrated for antiferromag- 
netism; it is not possible to have perfect antiferromagnetic order on it. In the 
afmll structure each atom has eight antiferromagnetically oriented and four 
ferromagnetically oriented neighbors, maximizing the antiferromagnetic inter- 
actions. In the case of a nearest-neighbor (n.n.) Heisenberg model, which has 
energy E = JiY.n.n.'^i ' the energy is independent of the angle between 
the moments of one antiferromagnetic pair and another (see Fig.4). If non- 
neigherest neighbor interactions are important, or if the Heisenberg model 
does not completely describe the energetics, the energy might be further low- 
ered if the antiferromagnetic pairs are oblique or perpendicular to each other 
(i.e. a ^ 0). 

The collinear afmll stru c ture h as 4 atoms per unit cell in space group Pmma 
( Steinle- Neumann et al. . 19991 ). In order to tile the lattice with the pattern 



shown in Fig.4 a cell with 8 atoms is obtained with space group A value 
of 1.6 was used for c/a for all of the hep based calculations here, which is 
close to the ground state value for the volumes studied here. Tests showed 
that c/a does not vary significantly (j.005) with magnetic state. Figure 5 
shows the energy versus angle a for different volumes. In all cases the energy 
decreases when a is varied from zero, indicating that the ground state is non- 
collinear. The most stable state is with a=90°, which minimizes the local 
ferromagnetic interactions. Table 1 shows the results of the TB model for hep 
Fe with ferromagnetic, admll, and the non-collinear structure with q;=90° 
(ncl). 

The ground s tate non-coll i near t otal energies of e -e were fit to a Vinet equa- 



tion of state (| Vinet et al.L |1987|) as a function of volume. Table 2 shows the 
results both using K'=4 and allowing K' to vary. For non-magnetic Fe, K'=4 
is not a bad approximation, but when magnetism is included K' >6. The 
results show that the TB model is in god agreement with the self-consistent 
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Fig. 2. Variation of magnetization energy versus wavevector for 7 — Fe for (0,0, q), 
with q in units of 27r/a , where a is the lattice constant). Wavevectors q = and 
q = 1 correspond to ferromagnetic and antiferromagnetic structures, respectively, 
q = 0.5 is a non-collinear spin state wit h relative angle betwe en the neighbor- 
ing spins of 7r/2. (a) LMTO resul t s from iMrvasov et all (|l992l ^. (b) LAPW re- 
sults from ISiostedt and NordstromI ()2002l l. (c) Results from this study using the 
best-fit value of / (1.017, 1.006, 1.003, 1.000, and' 0.993 eV for the volumes shown), 
(d) Results from this study using I =0.94 eV. The tight-binding results show 
the correct behavior of ferromagnetic at large volumes and antiferomagnetic at 
si nall volumes (high pressures) w ith a non-collinear transition region. The results 
of'Si5ste dt and Nordstr"5^ ()2002h seem more different than our results and those of 
livlrvaso v et al.l (Il992h than can be explained by the atomic moment approximation. 
Note that the line dashes and volumes correspond for (a), (c) and (d) but not for 
(b). ^ 




Fig. 3. Var iation of moments ver sus wavevector for 7 — Fe for (0,0,q). (a) LMTO 
result s from lMrvasov et aL ( 19921 ) . (b) LAPW results from Siostedt and NordstromI 
{ 2m± . (c) Results from this study using the best-fit value of / (1.017, 1.006, 1.003, 
1.000, and' 0.993 eV for the volumes shown), (d) Results from this study using 
/ =0.94 eV. Note that the line dashes and volumes correspond for (a), (c) and (d) but 
not for (b). A high-spin low-spin transition is evident in all cases for ferromagnetic 
and low-g spin- waves with increasing pressure. The tight-binding model apparently 
tends to give too small a moment in the low spin regime, but is in generally good 
agreement. 
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Fig. 4. The motif for an hep lattice with antiferromagnetic interactions. One pair 
of antiferromagnetic iron atoms is at an angle a to another pair. In the Heisenberg 
model with near-neighbor interactions, the total energy ins independent of a. The 
coUinear afmll structure is represented by a = 0. 
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Fig. 5. The motif for an hep lattice with antiferromagnetic interactions. One pair 
of antiferromagnetic iron atoms is at an angle a to another pair. In the Heisenberg 
model with near-neighbor interactions, the total energy ins independent of a. The 
collinear afmll structure is represented by a = 0. 

LAPW computations of ISteinle- Neumann et al. Including magnetism 

lowers the bulk modulus, and including non-coUinear magnetism lowers it fur- 
ther. The resulting equation o f state fo ncl is in very good agree r nent w ith the 
experiments (Fig. 6,Table 2) ilephcoat et all Il986l : Mao et al.L Il99nh . When 
magnetization is not included the disagreement between the experiment and 
theory is 75% and 9% for the bulk modulus and the equilibrium volume re- 
spectively. Including the afmll structure as the ground state of hep reduces 
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Table 1 

Magnetization energy (mRy/atom) for different atomic volumes V (bohr^). The 
value of Stoner parameter I (eV) for the different volumes is given in the first 
column. The magnetic moments are given in parenthesis. The non-collinear (ncl) 
structure is for a=90° 



V (au) 


I(eV) 


Ferro 


afmll 


ncl 


90 


1.01 


27.1 (2.6) 


17.7 (2.28) 


21.6 (2.38) 


80 


1.05 


7.2 (2.6) 


7.1 (1.79) 


9.7 (1.99) 


75 


1.1 


0(0) 


3.1 (1.47) 


4.6 (1.67) 


70 


1.14 


0(0) 


0.55 (0.93) 


1.1 (1.13) 


65 


1.19 


0(0) 


0.08 (0.33) 


0.02 (0.19) 




V (bohr') 



Fig. 6 . Equation of state of e-Fe. Symbols are experiments from I jephcoat et al 



19861 ) (open) and iMao et al.l (|l990|) (solid). The solid and the dashed curves are 
our magnetic and non-collinear NC 1 magnetic theoretical tight-binding results. 
Inset shows the magnetization energy per atom and moments in NC 1 versus 
volume. The upper non-linear scale shows the corresponding theoretical 

pressures. 

the disagreement to 25% and 5% for bulk modulus and equilibrium volume 
f Steinle- Neumann et al. . 19991 ). When the non-collinear spin state energies 
were used in the equation of state the bulk modulus was within 10% of the 
experiment. 



In spite of theoretical evidence for magnetism in e-Fe, and the great improve- 
ment in the equation of state when magnetism is included, the experimental 
situation i s unclear. Mo s sbauer spec troscopy shows no evidence of magnetism 
in hep Fe ( Tavlor et al. . 1982L 199l| l. X-ray absorption spectroscopy also has 
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Table 2 

Comparison of experimental and theoretical values of equilibrium volume (Vq) and 
bulk modulus (Kq) for e-Fe. 



be (GGAj 


Vq (bohr ) 




K 


bxDtdViao et al.. 1990) 


75.4 


Ibo 




JNon-JVlametic iGGAl I htemle-JNeumann et al.. 1999) 


69.0 


292 


A A 

4.4 


JNon-JVlagnetic (this study) 


69.1 


300 


4.0 


Non-Magnetic (this study) 


68.8 


297 


4.6 


Collinear (afmll) (Steinle-Neumann et al.. 1999) 


71.2 


209 


5.2 


Collinear (afmll) (this study) 


71.3 


240 


4.0 


Collinear (afmll) (this study) 


70.1 


214 


6.3 


Non-CoUinear (this study) 


72.0 


227 


4.0 


Non-CoUinear (this study) 


70.7 


195 


6.5 



been interpreted to show lack of magnetism (jR.ueff et al.L Il999l ) , but due to 
the absence of an absolute calibration, and sensitivity of the spectrum to 
changes in the density of states, the experiments show only that magnetic 
moments are lower in e-Fe than a-Fe, a result consistent with theory. On the 
other hand, there is some independe nt evidence of magn etism from Raman 
spectroscopy, which shows two peaks ( Merkel et ah . 200(1 ). instead of the one 
expected in the hep structure. The frequencies and splitting of these peaks 
is pre dicted well from first-principl es calculations for magnetic ordered afmll 
e-Fe f Steinle-Neumann et ah . 20031 ). and the absence of observed splitting in 
Mossbauer will be explained in a subsequent paper. 

Thanks to B. Fultz, A. Goncharov, R. Hemley, A. Jephcoat, L Mazin, H. 
Olijnyk, W. Pickett, and G. Steinle-Neumann, L. Stixrude, V. Struzhkin for 
helpful discussions. S. Mukherjee performed technical assistance. This work 
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W-7405-ENG-48 and the National Science Foundation EAR-998002. Com- 
putations were performed on the Cray SVl at the Geophysical Laboratory, 
supported by NSF grant EAR-9975753 and by the W. M. Keck Foundation. 
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